Four-Photon (In)Distinguishability Transition 



Malte C. Tichy, 1 Hyang-Tag Lim, 2 Young-Sik Ra, 2 Florian Mintert, 1,3 Yoon-Ho Kim, 2 and Andreas Buchleitner 1 

1 Physikahsch.es Institut, Albert-Ludwigs-Universitat, 
Hermann- Herder-Str. 3, D-79104 Freiburg, Germany 
2 Department of Physics, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, Korea 
s Freiburg Institute for Advanced Studies, Albert-Ludwigs-Universitat, Albertstr. 19, D-79104 Freiburg, Germany 

(Dated: March 15, 2011) 

We demonstrate the conspiration of many-particle interferences of different degree to determine 
the transmission of four photons of tunable indistinguishability through a four-port beam splitter 
array. The probability of certain output events depends non-monotonically on the degree of distin- 
guishability, due to distinct multi-particle interference contributions to the transmission signal. 

PACS numbers: 03.65.-w, 05.30.Jp, 05.30.-d 



The symmetrization postulate imposed by the indis- 
tinguishability of particles is a fundamental quantum 
concept with no classical counterpart, which strongly 
influences the behavior of matter at any energy scale. 
In quantum optics, the most prominent manifestation 
thereof is the Hong-Ou-Mandel (HOM) effect [T]: two 
indistinguishable photons falling on the opposite input 
ports of a balanced beam splitter always leave the setup 
together, i.e. all amplitudes with two photons in dif- 
ferent output modes interfere destructively, while events 
with both photons in the same mode are enhanced due 
to the photons' bosonic nature. When tuning the transi- 
tion from distinguishability to indistinguishability of the 
photons, which can be described by a single parameter 
[2], the visibility of the HOM dip in the probability to 
detect one photon per mode increases monotonically. 

Signatures for the full indistinguishability of more than 
two particles can be observed, e.g., when many indistin- 
guishable photons bunch at one output mode of a beam 
splitter [5H5]. In addition, when several particles enter a 
multiport beam splitter [5] simultaneously, many-particle 
interferences strongly influence the probability of indi- 
vidual counting events, such that many distinct events 
with a given number of particles per output port are to- 
tally suppressed [7HS]- The transition between fully dis- 
tinguishable and fully indistinguishable particles has not 
received attention beyond its impact on bosonic bunch- 
ing [4]. Its role for the behavior of other events, e.g. the 
ones that are suppressed for indistinguishable particles 
[HI IH] i is widely open. A thorough understanding of par- 
tial indistinguishability is, however, mandatory for the 
experimental characterization of the degree and nature 
of many-particle interferences. 

In the present Letter, we consider four photons that 
propagate through a four-port beam splitter array. In 
contrast to an intuitive extrapolation of the features of 
the well-known two-photon case pQ , and of the bunching 
behavior of many photons [I] , we show that (i) the degree 
of distinguishability manifests itself in non-monotonic 
event probabilities, and that (ii) events with large occu- 
pation numbers are not necessarily more likely with in- 



creasing indistinguishability, despite the bosonic nature 
of photons. These results are established by quantitative 
predictions on experimentally directly accessible quanti- 
ties. 

We study many-particle interference within a setup 
with four ports [10], see Fig. [TJa). In this arrangement, 
no single-particle (Mach-Zender-like) interference can oc- 
cur, what makes it perfect to study the genuine manifes- 
tation of many-particle interference. 




FIG. 1: (color online) (a) Diamond-shaped multiport (bal- 
anced) beam splitter array. Four ingoing modes a\ are redis- 
tributed onto the output modes St. The path length Xj of 
each incoming mode controls the mutual distinguishability of 
the particles in the setup, (b) Creation of four photons by 
spontaneous parametric down conversion, by double-passage 
of a laser pulse through a non-linear crystal. 

Single-particle evolution is here described by a unitary 
matrix which relates particle creation operators of input 
and output ports, a\ u and b\ U1 , respectively, via 
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with the photon frequency uj unchanged, and U(a, 4>) 
given by 

UM) = - \ \ \ \ . (2) 
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For any value of a and <f>, the matrix U(a, 4>) is a complex 
unitary matrix with \Ujk\ = 1/2, i.e. a Hadamard matrix 
[IT] . The phases a, have two rather distinct physical in- 
terpretations: <fi absorbs all relative phases between input 
modes, and is well-known from the two-mode interfer- 
ence of many particles, e.g. in NOON-state interferometry 
[12]. a is of distinct origin: It corresponds to the phase 
enclosed by the interfering modes (similar to a Sagnac 
interferometer [H]), emerges only in the four- mode case 
[IT] , and effectively controls the relative phases of the 
output components, conditioned on the input modes. 

Four photon states are created by the double passage of 
a laser pulse through a non-linear crystal, through spon- 
taneous parametric down conversion, see Fig.Qb). Such 
states consist of the coherent superposition of a quadru- 
plet part, with one photon in each mode, created when 
the laser pulse induces one pair of photons at each time 
it passes the crystal, with two double-twin parts, where 
four photons are distributed among two modes, result- 
ing from events where the pulse generates two pairs in 
the one or in the other passing direction TBI . Since all 
three processes occur with the same probability, the ini- 
tial state reads 

/ 4 ~f 2 -t 2 -T 2 -t 2 \ 

I*) = ^ IK, + + a ^f^\ |0) . (3) 

The operator 

creates a single photon with central frequency loq and 
spectral width Aw at the input port j, at time tj. The 
arrival times tj of the photons can be tuned through vari- 
able path lengths Xj = c • tj depicted in Fig. [IJa), and 
hence their overlap, or indistinguishability, by virtue of 
I (0| Gi.tXtfc l°> I 2 = exp(-±Aw 2 (^- - t k ) 2 ) A change of 
relative path lengths induces a phase-shift between the 
input-modes, which can be accounted for in <fi: 

4>^ 4> + ui(ti + t 2 -t 3 - u) (5) 

All possible final states that emerge from the multi- 
port can be characterized in terms of the photon number 
detected in each port. For four photons, simple combina- 
torics yields 35 distinct events which can be labelled with 
vectors s = (si, S2, S3, S4), where Sj is the number of par- 
ticles in port j. Their order is hereafter given by their rel- 
ative abundance in the fully distinguishable case, eq. ^ 
below, and such that vectors s with large si come first, 
i.e., Si = (4, 0,0,0), 3, = (0,4,0,0),...,s 5 = (3,1,0,0), 
and, finally, s 35 = (1, 1, 1, 1). 

For the interpretation of the interference effects that 
we will discuss hereafter, it is useful to decompose the 
initial state ([3l into a superposition of states in which 



the photons in different modes are all either fully indis- 
tinguishable or fully distinguishable. In the two-photon 
case, only two contributions (distinguishable and indis- 
tinguishable) arise. For four photons, the Gram-Schmidt- 
orthogonalization of the single-particle states that ap- 
pear in ([3]) yields 4! — 24 different terms for the four- 
particle state, since the description of the jth particle 
needs up to j non-vanishing components. Each term 
corresponds to a certain distinguishability setting that 
we will denote by {ii, 13, 14}, where photons in port 
k and I are indistinguishable if i/. = %i- Fully in- 
distinguishable/distinguishable particles correspond to 
{1, 1, 1, 1}/{1, 2, 3, 4}, and only indistinguishable parti- 
cles interfere. Hence, the distinguishability setting deter- 
mines the degree of many-particle interference. In gen- 
eral, a setup with given arrival times for all photons cor- 
responds to a situation in which several distinguishabil- 
ity settings contribute to the initial state, with different 
degrees of multi-particle interference occurring simulta- 
neously. 

The extreme case of full distinguishability {1, 2, 3, 4} is 
realized when all photons have pairwise delays — tj\ 3> 
1/Aw. In this case, we can safely neglect the expo- 
nentially suppressed components that still exhibit many- 
particle interference, and simple combinatorics can be 
applied to yield the output event probabilities (remem- 
ber that no single-particle interference can occur in our 
setting) : 

The opposite, fully indistinguishable limit {1,1,1,1} is 
realized for t\ = t 2 = t 3 = i 4 , when all photons can 
interfere perfectly. 

The resulting probabilities for five representative 
events are compared in Tab. [TJ for the fully distinguish- 
able and indistinguishable case, respectively. The event 
s*x4 = (0,1,0,3) is fully suppressed for indistinguishable 
photons, for any choice of the phases a, <j>, while $21, S22 
and S35 exhibit an intricate a, (^-dependence. Most im- 
portantly, there is no unambiguous correlation between 
the event probability and the photon distribution on the 
output modes: as anticipated above, si = (4, 0, 0, 0) may 
be enhanced (as expected for bosonic bunching) as well 
as strictly suppressed (for the experimentalist's choice 
6 = 7r). Unexpectedly, S35 = (1,1,1,1) may be en- 
hanced up to a weight 32/9. In stark contrast to the 
two-photon HOM effect, the manifestation of perfect in- 
distinguishability is not unique in the present, multi- 
particle interference scenario, and has no intuitive inter- 
pretation in terms of the occupation of modes. This is 
highlighted by Fig. [2j which displays all event probabil- 
ities, for fully distinguishable and indistinguishable par- 
ticles, and three different choices of a and 0, and thus 
demonstrates the loss of any indistinguishability-induced 
hierarchy in the event probabilities. Furthermore, not 
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Event 


fdist 


Pid(s) 


Si 


= (4,0,0,0) 


1/256 




cos 4 (0/2)/8 


Sl4 


= (0,1,0,3) 


1/64 







S21 


= (0,2,0,2) 


3/128 


[cos (a) 


+ cos(a + 4>)] 2 /48 


S22 


= (0,0,2,2) 


3/128 


[1- 


3cos(2a + 0)] 2 /48 


S35 


= (1,1,1,1) 


3/32 


[cos (a) 


- cos(a + 4>)f /12 




TABLE I: Event probabilities for fully distinguishable (Pdist) 
and indistinguishable particles (fid)- 



only are the event probabilities no unambiguous wit- 
nesses of (in)distinguishability any more, but they even 
evolve non-monotonically with decreasing distinguisha- 
bility of the particles, as can be demonstrated in our set- 
ting, by subsequently rendering different pairs of photons 
indistinguishable, and thus adding interference terms be- 
tween an increasing number of photons. 
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FIG. 2: (color online) Event probabilities of all 35 possible 
output configurations Sj , for distinguishable (black diamonds) 
and fully indistinguishable particles, with (f> — a = (blue 
squares), (j> = 0,a — 7r/4 (red circles) and <f> = n/4,a — 
(brown triangles). 

Let us focus on the effects of partial interference on 
the three events s 35 = (1, 1, 1, 1), S14 = (0, 1, 0, 3), S21 = 
(0,2,0,2) (all listed in Tab. [I]), within two exemplary 
distinguishability transitions: For specificity, we assume 
that the single photons have a central wavelength of 
Ao = 780 nm, and a full width at half maximum (FWHM) 
of Aa = 5 nm, what corresponds to a coherence length 
of l c fa 122 /xm. 

(I) We first adjust the path lengths in a continuous 
way, by parametrizing x% — 0, X2 — y, x% — —y and 
x± = 2y. Hence, for y = —180 /im, the particles are 
fully distinguishable {1, 2, 3, 4}, while for y = fim, their 
wave functions fully overlap, the particles are fully indis- 
tinguishable {1, 1, 1, 1}. During the transition, several 
distinct contributions to the output signal arise simul- 
taneously (see Fig. [3^,). We initially choose a = and 
= 0. The latter parameter remains constant for all 
values of y with our chosen parametrization (cf. (|5]))- 
Consequently, only variations on length-scales of the or- 
der of the coherence length l c appear. 



The signal evolution during the (controlled) indistin- 
guishability transition is shown on the bottom left of 
Fig. [3j The event s 35 = (1,1,1,1) is progressively sup- 
pressed as higher order interference contributions, from 
two- over three- to four-photons, kick in. This contrasts 
with the clearly non-monotonic event probability of si. 4 — 
(0, 1,0,3). While two- and three-photon-interference in- 
crease its probability, it is eventually strictly suppressed 
for fully indistinguishable particles, resulting in a maxi- 
mum for partial distinguishability. For s 2 i = (0,2,0,2), 
a step-like-behavior results from the competition of de- 
structive three-particle- and constructive four-particle in- 
terference. 

(II) Another remarkable impact of the distinguisha- 
bility transition is born out when we adjust the path 
lengths step-by-step: Starting from fully distinguishable 
photons, with x\ = 0, x-i = 220 /im, x$ = 440 /jm, 
X4 = 660 /im, we first tune down X2, and subsequently 
X3 and x i} as indicated in Fig. [3^c) , with a = n. 
The thus chosen, dominant distinguishability settings are 
{1,2,3,4} -> {1,1,3,4} -> {1,1,1,4} -> {1,1,1,1}, and 
this step-wise transition is reflected by the event prob- 
abilities in Fig. [3^d). In particular, the probability of 
s*x4 = (0, 1, 0, 3) is again non-monotonic, due to a destruc- 
tive two-photon contribution, followed by a constructive 
three-photon contribution for the {1, 1, 1,4} setting. 

An additional effect manifests when four-photon inter- 
ference sets in: the event probabilities start to depend 
on the phase 4>. This leads to fast oscillations on the 
scale of the photon wavelength A, an intrinsic feature of 
the interference of four or more photons, shown in the 
plots as shaded area between minima and maxima. The 
dependence of the event probability on <f> for fully indis- 
tinguishable photons (Tab. [I]) explains the onset of fast 
oscillations of the event probabilities for S21 and S35, and 
the growth of their amplitudes with increasing {1, 1, 1, 1} 
contribution. Unexpectedly, also si 4 = (0,1,0,3) ex- 
hibits such oscillations which, however, vanish when four- 
photon interference fully dominates, since this event is 
then strictly suppressed (Tab.[l]). It is the interplay of the 
{1,1,1,4} and {1,1,1,1} settings that leads to such de- 
pendence on </>: the double-twin part of the wave function 
in the setting {1,1,1,4} interferes with the quadruplet 
part in the {1,1,1,1} setting, resulting in four-photon 
interference that depends on (f>, and a local maximum of 
the event probability at the point where the contribu- 
tions of {1, 1, 1, 4} and {1,1,1,1} are equal. This feature 
is specific to the partial distinguishability of the photons. 

The discussed phenomena constitute show-case exam- 
ples for the intricate effects that manifest when quantify- 
ing or controlling the particles' degree of indistinguisha- 
bility in an experiment. Despite the mere doubling of the 
number of particles with respect to the HOM setup, the 
(in) distinguishability transition cannot be explained any 
more by extrapolation of the two-photon effect: interfer- 
ence dominates over bosonic bunching, and events with 
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FIG. 3: (color online) Evolution of event probabilities during the transition from distinguishable to indistinguishable particles 
(bottom row). Left: continuous transition, parametrized by y, with a = 0. Right: step-wise transition, with corresponding 
interferometer path lengths, Fig. [Tj on top of panel (c), with xi = /im, a = 7T. (a,c): Distinguishability type contributions 
to the output signal: While for the step-wise transition, at most two distinguishability settings are relevant at a time, several 
may contribute in the continuous transition. (b,d): Event probabilities of S35 = (1,1,1,1) (blue), S14 = (0,1,0,3) (red) and 
S21 = (0,2,0,2) (green). The shaded areas in (d) represent four-photon interference fringes oscillating as a function of (j>. 



large occupation numbers are not necessarily enhanced 
when approaching indistinguishability. The phase a 
which quantifies the phase enclosed by the setup con- 
stitutes an input- and output-mode dependent parame- 
ter with measurable impact, as a physical consequence 
emerging from the formal definition of four-dimensional 
complex Hadamard matrices [TT]. Finally, not only the 
degree, but also the quality of interference changes with 
the number of interfering particles. 

Our results impressively demonstrate that we need to 
fully abandon the idea that many-particle interference 
manifests itself in a unique fashion that can be predicted 
from the bosonic nature of particles alone, as usually 
stated in the two-particle case [15] . Constructive and 
destructive interference effects occur for all final events, 
and they can be turned into one another by variation of 
the phases. Our decomposition of the initial state into 
distinguishability settings does not only offer a powerful 
tool for the computation of event probabilities, necessary 
for the experimental implementation of many-particle in- 
terference, but also allows us to understand and clas- 
sify the occurring phenomena. Due to the combinatorial 
explosion of possible many-particle paths, a description 
of the indistinguishability transition becomes prohibitive 
for many more than four particles, and it remains to be 
studied whether one can find other, coarse-grained ob- 
servables [14] with a monotonic dependence on a single 
parameter, to quantify many-particle distinguishability. 
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